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We present the results of a study of specific heat on a single crystal of Pro.63Cao.37Mn03 performed 
over a temperature range 3K-300K in presence of and 8T magnetic fields. An estimate of the 
entropy and latent heat in a magnetic field at the first order charge ordering (CO) transition is 
presented. The total entropy change at the CO transition which is ~ 1.8 J/mol K at 0T, decreases 
to ~ 1.5 J/mol K in presence of 8T magnetic field. Our measurements enable us to estimate the latent 
heat Leo ~ 235 J/mol involved in the CO transition. Since the entropy of the ferromagnetic metallic 
(FMM) state is comparable to that of the charge-ordered insulating (COI) state, a subtle change in 
entropy stabilises either of these two states. Our low temperature specific heat measurements reveal 
that the linear term is absent in 0T and surprisingly not seen even in the metallic FMM state. 
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The discovery of a number of fascinating properties 
like colossal magneto- resistance, charge / orbital order- 
ing and electronic phase separation in manganites with 
generalized formula Rei^ x Ae x MnO^ (Re being a triva- 
lent rare-earth ion, Ae a divalent alkaline earth element) 
has resulted in a spurt of research activities For 
certain values of x, close to 0.5, these manganites undergo 
a first order transition at certain temperature Tco to a 
Charge Ordered Insulating (COI) state where the Mn 3+ 
and Mn 4+ species arrange themselves in a commensu- 
rate order in the lattice. The charge ordering transition 
in these oxides is accompanied by a large change in vol- 
ume and hysteresis in resistivity and is believed to be a 
first order transition. A fascinating aspect of the COI 
state (which is also accompanied by orbital ordering) is 
that it is unstable under an applied magnetic field and 
there is an insulator- metal transition (melting) of the 
COI state to a ferromagnetic metal state (FMM) below 
a temperature Tmh- 

In this paper we have investigated the specific heat and 
related thermodynamic quantities in a single crystalline 
CO system over a wide temperature range (3K < T < 
300K) and in a magnetic field upto 8T which can melt 
the CO order. A study of specific heat over an extensive 
temperature range is crucial in understanding the nature 
of the CO transition and in addition it can provide values 
of various fundamental parameters of manganites like the 
density of states at the Fermi level N(Ep), the Debye 
temperature #d, ferromagnetic / antiferromagnetic spin 
wave stiffness constant etc. 

Measurements of specific heat in manganites particu- 
larly with CMR composition has been reported before, 



particularly in the low temperature region. In the CO 
systems there are reports of specific heat measurements 
both in single and polycrystalline samples mainly at low 
temperatures. Specific heat measurements in CO sys- 
tems near Tco have been done in polycrystalline samples 
and no magnetic field data have been reported. There 
exists no data on thermodynamics associated with the 
transition at Tmh- 

The charge ordering transition in these oxides is be- 
lieved to be a first order transition, as mentioned before. 
Measurements across the CO transition should therefore 
show a discontinuous jump in entropy as the compound 
absorbs latent heat from the bath to transform to a new 
phase. However, till date there exists no proof from 
calorimetry that the transitions are indeed first order and 
whether a latent heat is released at the transition. There 
has been no calorimetric investigation of the melting of 
COI state to a FMM state, the transition which is also 
believed to be a first order transition. 

In this paper we investigate this fundamental issues 
related to changes in entropy and latent heat across the 
charge-ordering transition near Tco and Tmh in a single- 
crystal of Pr .63Cao.37Mn03. Pri-^ Ca^MnOa happens 
to be a prototypical and perhaps the most studied charge- 
ordered system. Due to its low tolerance factor it re- 
mains insulating for all values of x. For x = 0.37 com- 
position, charge and orbital ordering occurs at Tco (~ 
235K), while a long range AFM order sets in only below 
Tjv (~ 170K). A small ferromagnetic component appears 
with the canting of AFM spins at a lower temperature, 
T CA (« 30K). 

Though there have been a number of studies on low 
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temperature specific heat of this particular system (T < 
20K) , a clear picture is yet to be emerge on the presence 
(or absence) of a linear term in specific heat. Recently 
there have been reports of large linear contribution to 
specific heat and appearance of excess specific heat asso- 
ciated with charge ordering [|| in this system. Since such 
a linear term is often associated with a electronic contri- 
bution to the specific heat, appearance of this term in an 
insulating sample is intriguing. In this paper, we have 
also investigated the low temperature region and reached 
certain definite conclusions about the linear term. 

In our experiments we have specifically asked the fol- 
lowing questions: 

1. Is it possible to identify the relevant transitions in 
zero and finite magnetic fields through calorimetric 
measurements ? 

2. What latent heat is released and the entropy change 
across the first order transitions, one at Tco in zero 
magnetic field and the other at Tm h in 8T mag- 
netic field ? 

3. When the FMM state is obtained from the COI 
state in presence of a magnetic field, do we obtain 
a linear term (arising from electronic contribution) 
in the specific heat ? 

4. Do we see a linear term in specific heat in COI state 
as reported by some investigators ? 

The remainder of the paper is divided into two prin- 
cipal sections. In the first section after the experimental 
section, we present and discuss the experimental data in 
the region T > 50K which essentially encompasses the 
region where most of the transition occurs and this also 
happens to be the region where no experimental data 
have been reported in the past either in presence of mag- 
netic field or in a single crystal. The second major section 
refers to the data at low temperatures where issues like 
linear term in the heat capacity, the Debye term, spin 
wave contributions etc. are looked into. In this region 
there are past studies , as mentioned earlier and we com- 
pare our results on these material with results from other 
CO systems. 

I. EXPERIMENTAL 

The crystal used in our experiment has been grown by 
float zone technique and has been used in a number of 
previous experiments by our group The resistiv- 

ity (p) vs. T curve in zero field and in a field of H = 
8T are shown in fig. 1. The COI state can be melted 
to a FMM state by application of a magnetic field of 8T 
at Tmh ~ 90K. In presence of magnetic field the region 
Tmh < T < Tco is termed as "mixed charge order" 
(MCO) region, where the COI phase coexist with the 



FMM phase. We have used a semi-adiabatic heat pulse 
technique to measure the specific heat in a wide temper- 
ature range of 2K - 300K and in OT and 8T magnetic 
fields. 

II. SPECIFIC HEAT AT HIGH (T > 50K) 

In figure 2 we show the specific heat (C p ) data over 
an extended temperature range 3K - 300K for both 
and 8T. The most prominent feature is the sharp peak 
a t Tco both for H = and 8T. The peak is much larger 
and sharper compared to that seen in polycrystalline ma- 
terials [f| (which is also plotted in the same graph) . The 
peak at Tco retains its narrowness and sharpness at H 
= 8T while it shifts to lower T with dT co /dT « -1K/T. 
In figure 3 we show the region close to T/y. The graph 
shows the actual observed step like feature in C p (at 
H = OT) at Tjv- In the same graph we show the C exc 
after subtraction of the lattice contribution (the proce- 
dure for subtraction of the lattice contribution is given 
in the next sub-section). A peak in C' exc « 6J/moleK 
at Tn is visible. This feature is suppressed by appli- 
cation of a magnetic field of 8T. The entropy change, 
AS exc (T]\[), associated with this transition has been es- 
timated to be w 0.5 — 0.8 J /mole. K. Lee et. al. 
reported for Pro.6C'ao.4Mrt03(polycrystal) /S.S exc (T^) 

0.6J/mole K at T N =160K. Ramirez et.al. for 
iao.35Cao.65-^"03(polycrystal) obtained AS , ea;c (T/v) w 
U/moleK. It can be seen that all these are much less 
than that what one would expect from a complete spin 
ordering. The calorimetry data therefore points towards 
an incomplete spin ordering at Tjv- 

Another very interesting feature in our calorimetry 
data is a clear signature of the magnetic field induced 
transition at Tmh ~78K (see fig. 2 and 9). This small 
yet distinct feature, discussed in detail later on, is the 
first signature of field-induced melting in a calorimetry 
experiment. 

A. Estimation of the Lattice contribution 

In this subsection we describe the procedure to esti- 
mate the lattice contribution to specific heat, the back- 
ground on which the specific heat contribution by the 
other degrees of freedom add up. Proper estimation of 
the lattice contribution will thus allow us to get the con- 
tribution of CO and magnetic ordering to C p . We define 
Ce X c as . 

Ce X c — Cp Clattice (1) 

According to Dulong and Petit 's law, the limiting heat 
capacity at high temperature for a compound with r 
atoms per molecule is expected to be 3rR, where R is the 
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gas constant (juj). In our case (r — 5) this limiting value 
of the lattice (vibrational) heat capacity turns out to be « 
125 J / mole K. From our data we find that at T = 300K, 
observed C p is « 90% of the Dulong Petit value. For most 
published data on manganites, C p at room temperature 
reaches this value ~ 100 -120 J / mole K. This signifies 
that bulk of the contribution indeed comes from the vi- 
brational heat capacity in this temperature range. This 
contribution must be subtracted from the observed data 
to obtain the excess specific heat. 

We have obtained the lattice contribution Ci a tu C e by 
fitting C p in a region 40K < T < 150K using three stan- 
dard models as described below. The calculated Ci at u C e 
is then extrapolated to T > 150K. Since, we are in- 
terested in C p in the range 80K - 300K, the lower tem- 
perature for lattice contribution has been limited only 
to 40K. The upper limit of 150K was considered since 
it was lower than the AFM and CO transition tempera- 
tures. We have not estimated the background by includ- 
ing the data about Tco ■ We are of the opinion that there 
is structural transition/ modification associated with the 
CO transition . As a result it is not advisable to include 
the data above Tco m the estimation of lattice heat ca- 
pacity below Tco- 

The following models were used for estimating the 
background lattice contribution : 
The Einstein model @: 

C Einstein = 3rRJ2iai[xi 2 e x '/(e Xi - l) 2 ] ; where x 4 = 
Kvei 'KsTi. In this model, all the 3rN independent oscil- 
lators populate three optical modes in ratios a\ : a 2 : 03 
having Einstein frequencies hi/Et ■ The best fit to the data 
is given the three Einstein modes with hv/Ks = 145 K, 
410 K and 625 K. It is interesting to note that a recent 
Raman measurement on x = 0.37 composition observed 
optical modes at w 360 K, 417 K, and 648 K. The 417K 
and 648K modes nicely matches with those seen in our 
specific heat measurements. Jl4| 
The Debye model O: 

C Debye = 9rR/x 3 D £ D x%e x °/(e XD - l) 2 dx D ; where 
xjj = hvo/KsT. The specific heat is due to collective 
low-frequency oscillations of phonons with a cut off fre- 
quency given by Iivd/Kb- For our sample, the best fit 
was obtained with 9o = 470-KT. This is typically the value 
of 9d seen in most oxides. 
The Thirring model Jl|: 

Crhirring = 3ri? J2n°=o B nU~ n ; where u = [(T/T b ) 2 + 1] 
The harmonic portion of the lattice specific heat can 
be expressed in a series with the above form where 
Tb w 0d/27t p7| . The above expansion permits the 
harmonic portion of the lattice specific heat to be fit- 
ted reasonably well down to temperatures ~ 50if even 
when the Debye temperature is ~ 500K. In our case, 
Tb = 65K and we used n upto 50. The Tb gives an esti- 
mated 6d ~ 410.R', which is close to that obtained from 
fitting the data with Debye model. 



In figure 4 we have shown the deviation of the observed 
data from the fit (i.e. AC/C — (C b s — C ca i c )/C b s ) in 
order to ascertain the extent of uncertainties involved in 
the background subtraction. We find that for T < 10QK 
the Debye model shows large systematic deviation. The 
Thirring model shows a systematic deviation above 150K 
(not shown in graph) . We find that over the whole range 
the least uncertainty is shown by the Einstein model 
(maximum deviation from fit ±5%) and this is a random 
deviation. 

We have used the same CLattice for both H = and 
H = 8T. We found that if we fit the C p at H = 8T from 
50K to 140K to the above models as we have done for the 
C p at H = 0T case, we end up in getting essentially the 
same C 'Lattice using the Einstein model. There is a small 
(< 10%) systematic deviation (for the Einstein model) 
which is not numerically significant to affect our results. 
This is shown in the inset of figure 4. For compounds 
containing Pr 3+ Cp may contain a contribution coming 
from the crystal field. We have estimated this crystal 
field contribution or C x t a i the crystal field data available 
for iso-structural compounds PrNiO$, PrGaO^. We find 
that in the range of interest (100K < T < 300K) this 
crystal field contribution < 7% at 100 K and is < 2% at 
300K @. 

■B. Cexc 

In figure 5 we have plotted the value of C exc for both 
H = and H = 8T for T > 120K. The data show the 
C exc using all the three models. The C exc data in the re- 
gion close to Tco is truncated in figure 5, because in this 
region the C exc is very large. In this temperature range 
the differences in the C exc obtained after the subtraction 
of the lattice contribution is well within 1-2 J/mole 
K for the Debye and the Einstein model. The Thirring 
model shows systematic deviation for T > 150/-T . Con- 
sidering all the models we find that near 2jv the uncer- 
tainty is the largest since C exc is low and this can be as 
large as ±25%. However close to the CO transition when 
C exc shots upto sa 220 J /mole. K. The uncertainty in the 
background estimation jails below ±2%. The C exc data 
also show clearly that barring the small region around 
Tjv (where C exc ~ 10% of C p ) and near TcO (where 
C exc > 50% of C p ) the contribution by other degrees of 
freedom compared to the lattice contribution is negligible 
in the temperature range T > 50K. 

C. Entropy change near Tco in H = and H = 8T 

A phase transition is signaled by a singularity in a 
thermodynamic potential such as free energy If there is 
a finite discontinuity in one or more of the first deriva- 
tives of the free energy, the transition is first order. At 



3 



a first order transition one expects a discontinuous jump 
in the entropy. At the charge ordering transition, (which 
is conjectured as first order transition based on observed 
hysteresis in transport data on field and temperature cy- 
cling), our system is expected to absorb the latent heat of 
transformation to transform from charge-ordered phase 
to a charge-disordered phase as it is heated through Tco- 
The temperature should remain constant till this process 
is complete and the entropy change is given as : 

AS 21 = S 2 - Si = L 2 i/T (2) 

where S\,S2 are the entropies of phases 1 and 2, L 2 \ is 
the latent heat associated with the phase transformation. 
At a first order transition, since the entropy changes dis- 
continuously, specific heat is actually undefined. Because 
at the transition point, the temperature does not change 
when heat is applied one would expect a 6 function like 
behaviour. In reality, however, one rarely sees such a 5 
function like behaviour of C p because the transition can 
get broadened either by the process of measurement or 
by the quality of the sample. This makes a part of the ex- 
pected entropy change to become continuous and it shows 
up as a finite measurable specific heat. Interestingly, it is 
only very recently the expected features in specific heat 
in a first order transition has been seen experimentally in 
rare-earth compounds p2] . In the following discussion, 
we have to keep in mind the above special features since 
we are measuring specific heat near a transition which is 
expected to be first order. 

In figure 6, C exc /T is plotted as a function of T within 
a small interval of ± 25K around Tco for both H = and 
8T. In figure 7, we have shown the excess entropy S exc (T) 
associated with this transition. The entropy, S exc (T), 
for T > 200K is calculated by numerically integrating 
C exc /T as S exc (T) = j 2Q0 C exc /TdT. The lower limit of 
the integral is purely a matter of convenience. It is far 
removed from T/v and Tco and at T = 200K, C exc » 3% 
of C p and is negligible. 

We estimate the total change in entropy ASt by us- 
ing linear extrapolation of the entropy values from above 
and below to Tco as depicted in the fig. 7. ASt ~ 
1.8J/mol K at OT and it decreases to 1.5 J /mol K at 8T. 
This total change in entropy takes place over an interval 
of w 10 -15K around T CO - The magnitude of AS* T (0T) 
agrees well with that found in polycrystalline samples 
by Lees et. al. || by integrating the area under the 
peak in C exc /T. In polycrystalline Lao. 3sCao. e^MnOs 
Ramirez et. al [Q obtained ASt~ 5J /moleK. In poly- 
crystalline Yb.5Cao.5-/WnO3 our group had observed ASt 
rj 2.5 J /moleK ||. The total entropy change at the 
CO transition thus seems to be a fraction of what one 
would expect from an order-disorder transition. The er- 
ror in the estimation of the entropy change arising from 
the uncertainty in estimation of background Chattice is 
not more than 5%. In a narrow temperature range over 



which the transition occurs the change in entropy being 
the difference of two quantities with similar background, 
the error arising from the estimation of the background 
is not severe. 

The charge-ordered phase is expected to have a lower 
entropy than the charge-disordered insulating phase at 
T > Tco- This implies that the sample absorbs latent 
heat to transform from the charge-ordered phase to the 
disordered phase. We would expect that atleast a part 
of the entropy change ASt is released as latent heat is 
absorbed from the total heat AQ supplied during the 
heat-pulse experiment and the sample temperature re- 
mains quite constant throughout the process. Since the 
change in sample temperature (AT) is quite small, the 
specific heat C p = (dQ/dT) p shows a very sharp peak 
near Tco- As discussed before at a first-order transition 
AT — > and ideally C p should be a 6 function at the 
transition temperature. 

The peak however is broadened due to two reasons : 

(a) Crystal quality : In most cases the expected S func- 
tion gets broadened by sample quality. The real crys- 
tal contains some defects / inhomogeneities which would 
lead to broadening of the peak. With improved qual- 
ity of the sample, the peak would be larger and sharper. 
For oxides containing multiple chemical constituents it is 
quite likely that such is the case. As shown in figure 1, 
where we have compared the C p of a ceramic pellet and 
a single crystal, the single crystal has a much higher and 
narrower C p at Tco as compared to the ceramic sample. 

(b) Measurement broadening : The second reason for 
broadening of the peak is connected to the measurement 
procedure itself. C p measurement involves measuring a 
finite temperature jump AT following the heat pulse AQ. 
Since C p is well behaved in the region away from Tco, h 
doesnot dependent much on the size of the temperature 
rise AT (as long as C p is not a very steep function of T). 
However, close to Tco, the height as well as the width 
ST of C p sensitively depends on AT. 

We have measured the latent heat by measuring C p 
with different AT as suggested by Gschneidner et. al. 
p2| . A sample of the data is shown in figure 8. We 
find that while away from Tco there is no dependence 
of C p on the rise of AT, as expected, the peak becomes 
narrower and higher when size of AT decreases as we 
approach Tco- This continues till AT < 0.5K. For this 
value we reach a limiting width (at half maxima) ST ~ 
2K. We believe that for AT < 0.5K, the width of the 
peak is not determined by the measurement but by the 
crystal quality. This observation sets a limit to the height 
and width of C p . We find from our measurement that the 
latent heat Leo /Tco — ASco ~ 1 J/mole K. Similar 
value of ASco was obtained for H = 8T. Our estimate 
of latent is a lower bound of the true latent heat and 
with a better crystal preferably with a transition width 
< IK, ASco I ASt will increase substantially and may 
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even — > 1, as expected for a strong first order transition. 
We note that such a narrow peak might arise from a very 
narrow second order transition that has been broadened 
by defects or inhomogeneities. However we find that the 
applied field shifts the position of the peak but doesnot 
appreciably broaden the peak. This we take as a proof 
of a first order transition at Tqo |2^1 ■ 

D. Entropy change near Tmh in 8T 

In a field of 8T the low temperature phase is an FMM 
or a spin-aligned metal, allowing for spin canting. But 
the phase is metallic as can be seen from figure 1. On 
heating, the FMM state (8T) becomes unstable towards 
the COI state formation and makes transition to the 
mixed charge ordered (MCO) state at Tmh- This is ac- 
companied by a jump in p. We would like to ask if there 
is a change in entropy associated with the melting of the 
COI state to FMM state at Tmh ? This question is inter- 
esting because both the phases have ordering of a kind, 
in addition the FMM phase is expected to have extra en- 
tropy due to presence of free electrons. A close look at 
the region T w 95K, as shown in figure 9, shows a small 
dip in C p which as explained below can be considered as 
signature of this melting. It is important to note that, as 
expected in zero field, this feature is absent. As shown 
in fig. 9, C P (8T) starts to show change at T = Tmh ~ 
88K where p shows a jump on heating and the specific 
heat transition is complete at lOOif after a small dip in 
C p , where the resistivity transition also stops. 

The small dip in C p is associated with a small heat 
release of m lOJ/mole at around 95K - 100K on heating. 
This suggests that on heating the FMM phase, stable in 
8T for low T, starts to disorder. This probably destabi- 
lizes the FMM phase with regions of high resistivity COI 
regions appearing in it. Eventually at T « Tmh the in- 
sulating regions increase in size and p shows a jump in p. 
The resistivity transition has an element of percolation 
associated with it and thus occurs at a different volume 
fraction of the new phase in contrast to the specific heat 
transition which occurs mainly when the bulk of the sam- 
ple is transformed. Nevertheless, the small but distinct 
change in C p close to Tmh is clearly seen. 

In the temperature range Tmh ^ T <C T^Qin a mag- 
netic field(i.e, the MCO region) we find an interesting ef- 
fect. The thermal relaxation time of the sample increases 
by more than one order of magnitude. This anomalously 
large relaxation is not observed at any temperature range 
in zero field or for T > T co and T < T MH in H = 8T. 
The existence of this large thermal relaxation necessi- 
tated that the sample be properly equilibrated before the 
data are taken. We have done that and the data pre- 
sented here are taken after the sample has been properly 
equilibrated thermally. The thermal relaxation, however, 



is interesting in its own right and has been discussed in 
a separate publication |25| . 

III. LOW TEMPERATURE C P IN AND 8T 
MAGNETIC FIELD 

As it has been pointed out before, that the issue of 
low temperature specific heat (T < 1QK) is controversial 
in manganites showing charge ordering, particularly in 
Pri_2;Ca x Mn03 system. A large number of observations 
on low temperature specific heat is available for differ- 
ent compositions of Vy\- x Ca^MnOs, mostly in poly- 
crystalline pellets and some in single crystals. 

In figure 10 we show the specific heat data of our sam- 
ple with that of different compositions as obtained by 
different groups @||. Such a comparison is meaningful 
because it brings out the essential similarity and differ- 
ences in these materials. A large linear term (7 ss 30.6 
mJ/mole-K 2 for x = 0.3 and 7 w 15.7 mJ/mole-K 2 for 
x = 0.35 ||) appears in polycrystalline samples with x 
= 0.3 and x = 0.35 composition which decreases as x — > 
0.5. We note that such a large linear 7 has not been seen 
any of the single crystal data. From a comparison of the 
data shown in figure 10 we can reach the following con- 
clusions: 

(i) The polycrystalline ceramic samples have higher spe- 
cific heat than the single crystal samples, 

(ii) There is a clear trend that the specific heat of the 
ceramic samples decreases as we approach x = 0.5. In 
particular, linear term in specific heat, is only observed 
for the ceramic samples and it decreases with increasing 
x. 

(iii) The specific heat of the single crystal samples with 
compositions x = 0.37 and 0.5 are very similar for 
T < WK. 

In figure 11 we have compared the specific heat of two 
close compositions: a polycrystalline x = 0.35 sample 
H and a single crystal sample with x = 0.37 composi- 
tion. The comparison is to elucidate how much of the 
excess specific heat is due to polycrystallinity and can 
we infer a likely origin of the excess specific heat. There 
are two mechanisms that can contribute to the excess 
specific heat: (1) Due to grains of small dimensions p7|j 
and (2)Due to two-level-systems(TLS) arising due to dis- 
order as in amorphous solids and several disordered crys- 
tals ^7j. If the excess specific heat, is due to grains then 
SC can be expressed as SC = C\T + C2T 2 , where the 
coefficients C\ and C2 are related to the average grain 
diameter < R > as < R >= ks [6fiv soun dCi}~ - 5 , where 
v sound is the sound velocity in the crystal p7| . We 
can fit the experimentally observed C po i y to the equa- 
tion Cp iy = C sc + SC. Using C sc from our observed 
data and the above expression of SC we obtain an es- 
timate of the average grain size < R > and we obtain 
< R >« 10~ 4 /im which is far too small compared to the 
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typical < R >> Ifim seen in most polycrystalline sam- 
ples. We can thus rule out finite grain size as the source 
of excess specific heat. 

If the excess specific heat 5C arises from TLS then it 
can be fit to a form SC = C\T + C 2 T 3 . This would arise 
from TLS with density of state P{E) = a + bE 2 , where 
a and b are constants j27j. Using this expression for 
SC we then fitted the observed C po i y as had been done 
before. The fit is shown in figure 11. From the fitted 
parameters C\ and C 2 , we arrive at values of a w 2.3 x 
10 35 erg" 1 / 00 an d 6 ~ 1.5 x 10 35 erg~ 3 /cc. These val- 
ues are quite comparable to but some what larger than 
those obtained in glasses. We conclude that extra specific 
heat of the polycrystalline samples arise from excitations 
which behave as TLS. We have no clear understanding of 
the origin of these TLS but we can speculate that they 
will arise from incomplete orbital /charge order that can 
happen in polycrystalline samples due to random strains. 
What is needed for TLS to occur is to have multitude 
of ground states with almost degenerate energy. This 
can actually be due to the predominant incommensurate 
charge ordering as seen in both x = 0.35 and x = 0.37 
compositions. As the charge-ordering becomes commen- 
surate with the lattice with x — > 0.5, the difference in 
the specific heat between the polycrystalline and the sin- 
gle crystal samples decrease. We thus suggest that the 
disorder in polycrystalline samples give rise to TLS low 
energy excitations typically with energy < 10K and as 
x -> 0.5, the P(E) -> 0. 

Next we attend to the low temperature data taken on 
our single crystal sample. Figure 12 shows the specific 
heat data for the x = 0.37 sample plotted as C p /T vs. 
T 2 , a customary way to plot the data in anticipation of 
a linear term. It is clear that the linear term is absent 
in the COI state. This is not surprising for an insulating 
sample, the linear term was also absent in the specific 
heat data observed in the ceramic samples of x = 0.4 by 
Lees et. al [Q. Neglecting the linear term, we have fitted 
the our observations to the following relation : 

C p = aT~ 2 + (3 3 T 3 + (3 5 T 5 (3) 

where, aT~ 2 hyperfine contribution caused by the local 
magnetic field at the Mn nucleus due to electrons in un- 
filled shells, /3 3 T 3 and P5T 5 are the lattice contribution to 
the specific heat, arising from phonons. A part of the T 3 
contribution to C p can also arise from AFM spin waves. 
The T 3 contribution is likely to be enhanced over and 
above the actual Debye contribution because of the pres- 
ence of AFM spin waves, since our COI sample has canted 
AF order in this temperature range. The results of fit- 
ting the data to the above equation are shown in table 1. 
The parameters obtained by us are very similar to that 
seen by Lees et. al 0. In table 1 we have collected the 
parameters from different published data for comparison. 

However, the interesting question is whether in a mag- 
netic field of 8T where the COI is melted into a FMM 



phase, we see appearance of a linear specific heat term. 
In a past investigation on polycrystalline samples of Lao. 5 
Cao.sMnOs in a magnetic field of 8.5 T |10| (which is 
high enough to melt the COI state) no linear term had 
been observed. The issue of absence of a linear term in 
the FMM phase obtained after melting the COI state 
is thus real. Particularly, when one compares it with 
the low temperature specific heat of the FMM phase as 
seen in the CMR region 7 w 7.8 mJ/mole K 2 in Lao.s 
Ca 2 Mn0 3 0, « 5.2 mJ/molc K 2 in La 7 Ca 3 Mn0 3 
& 

The Cp/T vs. T 2 for our sample in presence of an 8T 
magnetic field are shown in figure 12. The specific heat 
data have been fitted using the eqn. 1, with the param- 
eters as shown in table 1. Interestingly, the specific heat 
in presence of magnetic field is distinctly lower than the 
specific in 0T for all T < 4QK. Lowering of specific heat 
in presence of magnetic field has been reported Smolyani- 
nova et. al §. Since resistivity measurements indicate 
the sample to be in a metallic state, one would expect 
an extra linear term appearing in the eqn.l in presence 
of magnetic field. As seen in fig. 12, we couldnot detect 
the presence of any linear term even in presence of 8T. 
From the uncertainty in the data we find an upper limit 
of 7 < O.lmJ/moleK 2 , which is too small. 

It is apparent from table 1, that contribution of the 
/3 3 term is halved in presence of magnetic field. This 
appears to be the main reason for lower C p in H = 8T 
compared to that in absence of magnetic field. This can 
be explained as follows. The T 3 term contains an addi- 
tional antiferromagnetic spin wave contribution for the 
COI sample. For an antiferromagnetic spin wave spec- 
trum where E = Dq, the magnetic contribution to spe- 
cific heat is cx T 3 . The melting of COI state to a FMM 
state leads to collapse of this AFM order and thus there 
will be no AFM spin wave contribution to the specific 
heat. This would result in a decrease of the T 3 term. 
However, presence of FM order should give a ferromag- 
netic spin wave contribution cx T 3 / 2 . We have attempted 
fitting a T 3 / 2 term to the 8T specific heat data. However, 
we could not detect any J" 3 / 2 term. Thus, it appears that 
the specific heat as observed in presence of magnetic field 
doesnot have a ferromagnetic spin wave contribution. 

Briefly, we find that the specific heat at low temper- 
ature needs more investigation in single crystalline CO 
systems and in high magnetic field. The existing data 
(which includes too few single crystal data) both in zero 
field and magnetic field (high enough to melt the CO 
state) docs not allow us to reach much definite conclu- 
sion other than the fact that in the FMM phase obtained 
by melting COI state in a high magnetic field does not 
have a large enough 7 that is comparable to FMM phases 
seen in manganites showing CMR behavior. 
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IV. CONCLUSION 

Our study of specific heat in a CO system done over 
an extended temperature range of 3K - 300K in presence 
of and 8T magnetic field gives us a number of useful 
and new information regarding the thermodynamics of 
CO transition. This paper gives the first clear measure- 
ment of the entropy in a single crystal and also the latent 
heat (at Tco) m zero field and in a magnetic field. The 
latent heat released proves the strong first order nature 
of the transition, which had never been shown to be so 
in actual measurement of specific heat. An estimate of 
the numbers like the latent heat in and 8T field sets 
clearly the bound and scale within which any theoreti- 
cal models must work. It is a very important number 
for any phase transition which gives the CO transition 
(both in OT and 8T) a thermodynamic basis. Important 
is also the observation of a finite entropy change at the 
CO melting in a magnetic field. The result shows that 
a small but finite and qualitative difference in entropy 
exists between the FMM phase and the CI phase at the 
melting transition in the magnetic field where the former 
(FMM) phase (which is the low T phase) has a lower en- 
tropy than the other phase (which is the high T phase). 
Interesting also is the smallness of the entropy. The low 
temperature C p (T < 15K), show no linear dependence 
on T as expected from an insulating sample. Surpris- 
ingly, this electronic term is absent even when the COI 
state is melted to a FMM state by a 8T magnetic field. 
An attempt had been made to understand the observed 
discrepancy between the specific heat of a single crys- 
tal and a polycrystalline sample in the low temperature 
regime. 
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Figure captions : 

figure 1: The temperature dependence of resistivity 
of Pro. 63 Cao.37Mn03 in presence of and 8T magnetic 
fields. 

figure 2: C p of single crystalline Pri-^.Ca^MnOs (x 
= 0.37) over 3K- 300K in presence of and 8T magnetic 
fields. Also plotted the C p of the polycrystalline x = 0.4 
composition. 

figure 3: C p and C exc close to TV in H = 0T. 

figure 4: Deviation of C p from the fit for different 
models. Inset shows the same near Tmh in H = 8T. 

figure 5: C exc obtained after subtracting the lattice 
background (as determined by the three models) for both 
H = and 8T. 

figure 6: C exc /T around T co in H = and 8T. 

figure 7: Excess entropy S exc as calculated from nu- 
merically integrating C exc /T around Tco in H = and 
8T. 

figure 8: Varying sharpness and width of C p near Tco 
for two different temperature rise AT = 0.5 and IK. 

figure 9: C p and resistivity near the Tmh- The arrow 
marks the melting of the COI state. 

figure 10: A comparison of C p of Pri_ K Ca K Mn03 
in the low temperature regime as observed by different 
groups, x = 0.5(sc),0.45(ce),0.35(ce) @, x = 0.4(ce) g 
and x = 0.37(sc) is our sample. 

figure 11: The excess specific heat, SC, as seen in 
the polycrystalline sample arising due to presence of two 
level states. 

figure 12: The low temperature C p of single crys- 
talline Pr .63Cao.37Mn0 3 in H = and 8T. The solid 
line is the fitted equation. 

TABLE I. The fitting results for the C p of Pti- x C& x 
Mn03 as obtained by various groups. The units of differ- 
ent quantities are : a (mj K/mole), 7 (mj/mole K 2 ), (3% 
(mJ/mole K 4 ) and #5 (x 1CT 4 mj/mole K 6 ) 



X 

(Pri-zCa^MnOa) 


a 


7 


P3 


ft 


0.3 a 


63.0 


30.6 


0.30 




0.35 a 


56.0 


15.7 


0.39 




0.45 a 


28.0 


3.1 


0.31 




0.5 a 


22.0 


2.4 


0.26 




0.4 b 


28.0 




0.54 


4.5 


0.37 (H = 0T) C 


85.0 




0.54 


4.5 


0.37 (H = 8T) 


400.0 




0.16 


18.0 



a from Ref. 1 1 
b from Ref. [| 
c our sample 
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